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Abstract 

We present a new calculation of the renormalized HQET Lagrangian 
at order 0(l/mjo) and discuss the consequences of the BRST invari- 
ance of QCD and the reparameterization invariance of HQET. Our 
result corrects earlier, conflicting calculations and sets the stage for 
the calculation of the renormalized currents at order 0(1/ tuq). 
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1 Introduction 



Heavy Quark Effective Field Theory (HQET) |l| has been established as the 
theoretical tool of choice for the description of mesons and baryons containing 
heavy quarks || . This derives from the fact that it is a systematic expansion 
in inverse powers of the heavy quark with well defined and calculable coef- 
ficients. Furthermore, its realization of the spin and flavor symmetry of the 
low energy theory is a phenomenologically powerful tool. 

However, since the expansion parameter Ah/ttiq is about 0.12 (using a 
hadronic scale Ah m 600MeV and niQ = nib) leading order calculations are 
not sufficient for precision calculations. Terms of order 0(l/mn) including 
leading QCD corrections have to be under control. 

An indispensable prolegomenon to the calculation of renormalized matrix 
elements of currents is the renormalization of the Lagrangian. Unfortunately, 
two calculations with conflicting results || §] have been reported. In this note 
we present the result of a new calculation of the renormalized Lagrangian, 
which differs from the previous two. We will demonstrate that our result 
satisfies important consistency conditions that are violated by the earlier 
calculations. 

This note is organized as follows: in section |2] we introduce our operator 
basis. Our result for the anomalous dimensions is presented in section |3|. 
In section |] we will discuss the consistency of this result and compare it in 
section [| with earlier calculations. Finally, we present the renormalization 
group flow in section || and conclude in section \t\. 

The renormalized currents will be presented in a subsequent note [H and 
phenomenological applications will appear later 0. A more detailed discus- 
sion of technical matters will be presented in J7| . 

2 Operator basis 

The Lagrangian of HQET is defined by a systematic expansion of QCD in 
inverse powers of the heavy quark mass 

£ HQET = K, {M) K + ^ £ W + £ W + O(^) . 

(1) 

At order 0(l/mg) there is only one operator h v (ivD)h v , which is indepen- 
dent of the spin and flavor of the quark, resulting in the celebrated spin-flavor 
symmetry of HQET. 
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At order 0(1/ tuq) there are three independent operators. We shall use 
the conventional basis 



Of = h v (iD) 2 h v (2a) 
Of = 9 -h v a» x F, x h v (2b) 

Of =h v {ivDfh v . (2c) 

Below, the operators and Of will also be called kinetic and chromo- 
magnetic respectively. The operator of vanishes by the equations of motion 
(EOM) 

ivDh v \Q) = (3) 

for heavy quark states \Q). For the renormalization of the Lagrangian, it is 
not necessary to include this operator, if the EOM are used consistently. It is 
however needed as a counterterm in the renormalization of the heavy quark 
currents || and will be included here. In addition, inclusion of the operators 
vanishing by the EOM allows to make the relations following from reparam- 
eterization invariance explicit. Finally, the extraction of coefficients with the 
aid of symbolic manipulation programs is more straightforward in the full 
basis, while the work induced by the additional operators is insignificant. 

At order 0(1/ mq) there are thirteen independent operators. They are 
grouped in four classes. Two of the local operators do not vanish by the 
EOM. We will denote them collectively by O^ and choose them as 



Or = h v iD^(ivD)iD^h v (4a) 

Of = h v ia» x iD^(ivD)iD x h v . (4b) 

The five remaining local operators vanish by the EOM. We denote them by 
Ceom an d choose the basis 

Of = h v (ivD)(iD) 2 h v (4c) 

Of = h v {iDf{ivD)h v (4d) 

Of = h v (ivDfh v (4e) 

Of = - 9 -hSvD)a^F^ x h v (4f) 

Of = - 9 -h v o^F, x (ivD)h v . (4g) 
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In addition to the local operators, there are the time-ordered products of the 
lower dimensional operators. There are three of them that do not vanish by 
the EOM. They will be denoted f {2) . 

= ~T { [h v (iD) 2 h v ] [h v (iD) 2 h v ] } (4h) 
T}? = '{T { [h v (zD) 2 h v ] [Ka^F^K] } (4i) 
4 2) = *4~T { [h v a» x F^ x h v ] [Ka^F^K] } (4j) 

o 

Below, the operators and will also be called double-kinetic and 
double- chromo-magnetic respectively. Finally there are three more time- 
ordered products 7^ OM that vanish by the EOM 

r x ( 3 2) = iT { [h v (iD) 2 h v ] [h v (ivD) 2 h v ] } (4k) 
4 2) = f T { [Ka^F^K] [h v (ivD) 2 h v ] } (41) 
4 2) = \T { [h v (ivD) 2 h v ] [h v (ivD) 2 h v ] } . (4m) 

Below we shall refer to the operators vanishing by the equations of motion 
as EOM operators, for short. 

For calculational convenience, we have not chosen a manifestly hermitian 
basis, but the results presented below show that indeed only hermitian linear 
combinations show up in counter terms. 



3 Anomalous dimensions 

The most convenient approach to the calculation of anomalous dimensions 
uses the background field gauge || with gauge fixing term 

-1 {D»{V)A»f , (5) 

because only the renormalization constants of those operators that are man- 
ifestly invariant under gauge transformations of the background field have to 
be calculated. Therefore, only the divergent three-point functions have to be 
calculated to derive the anomalous dimensions. Furthermore, the Ward iden- 
tities for the classical background fields are particularly simple and provide 
a powerful tool for checking our results. 
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The background field gauge with an arbitrary gauge parameter allows 
an independent test of the consistency of our results by comparing the £ 
dependence with general results M for the one-loop effective action. Below 
we shall write £ = 1 — £ for the gauge parameter. 

The anomalous dimensions at order Oil/mq) are well known JL0J (the 
gauge parameter dependence has been calculated in |Tl||): 
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(i) 



(° 



v° 





\C A 




2d\ 



(6) 



where we have introduced the shorthand C\ = Cf(1 + £/2). Here and below, 
we have extracted the common loop factor a/n from (|J). 

The matrix of anomalous dimensions can naturally be written in block 
form, separating local operators from non-local operators and separating 
EOM operators from the rest: 
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(2) 



o(2) o 



(2) 

EOM 



fm f 

J- J-T, 



(2) 

EOM 



0(2) 




* (2) EOM 

Till 




] 


d (2) 

EOM 





* (2) EOM 
^1,2 




fm 




a (2) EOM 
7m, 1 


Jn 


. (2) EOM 
7nj 


<f(2) 
J EOM 


\ 


- (2) EOM 
7 m ,2 





. (2) EOM 
7n,2 / 



(7) 



The upper right block has to vanish, because Weinberg's theorem |T2[ guar- 
antees that the renormalization of the local operators does not require coun- 
terterms from the time-ordered products. The other three blocks display the 
block triangular structure required by the fact that the renormalization of 
EOM operators only induces counterterms that are EOM operators them- 
selves. 

The anomalous dimensions have been calculated manually. These cal- 
culations have been verified with the help of FORM [13| as a warm-up for 
the renormalization of the currents [[J, which requires the use of symbolic 
manipulation programs for economical reasons. 

We start the presentation of the results with the renormalization of the 
local operators 



7i 



(2) 



-\c A 

0, 



(8a) 
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~ (2) EOM 

7ii 



\C A \C A -2C F (l + 



fo -CV(l + 2£) 

-C F (l+2£) 



-\Ca 



~\C A , 



- (2) EOM 
7l,2 















V 



-CU o 
-C A J 



f8c) 



The renormalization of the time-ordered products requires local counterterms 
as well. The corresponding anomalous dimensions are 



(2) 
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a (2) EOM 
7m.l 



6 Ca — 3CV 



\ 



V -I^a y 







(8d) 



AF AF 





5 /m 5 /-» o/^' 



VT2 



— 2C F C A — 2C| 



\Ca 



using the shorthand C AF = j^Ca + CV(y + £) 



(2) EOM 
7m,2 



/^-C| -C* 6C F (l + 


-a 



-|C A + C| -±C A + C« 



(8f) 



The blocks ^fn \ ^ 2 j EOM and 7^ £jW ™ are given by the sum of the appropriate 
anomalous dimensions from 7W. 



, (2) EOM 



4 Consistency of the results 

The symmetries of HQET entail relations among the anomalous dimensions 
that can be used to check the result (§). Such consistency checks are useful 
in the present case and are of vital importance in the considerably more 
involved renormalization of the HQET currents ||. 
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4.1 Gauge invar iance 

As alluded to in the previous section, we have used the Ward identities for 
the background fields and the ^-dependence of the one-loop effective action 
for verifying our results. 

From the simple QED-like Ward identity for the background field 

q a T a (p,q) = S(p + q)-S(p), (9) 

with S and T a denoting the one-particle-irreducible two- and three-point 
functions with one operator insertion respectively, follows that the coun- 
terterms proportional to Ca (i.e. the non-abelian contributions) have to be 
transversal. Our result (§) passes this consistency check. The technical de- 
tails will be presented in J7|. 

In background field gauge, the ^-dependence of the renormalized effective 
action at one loop order is known to have the following form || 

-2P—f a = f (0) * f (1) 4- f (0) * F (1) 4- F (0) * f (1) + 
^ Q£ Aa Q JLh v h v ' h v A a h v MLh v h v A a h v MLh v ' " 

(10) 

where denotes the effective action at n-loop order, the same effective 
action with the gauge fixing term (||) subtracted and the effective action 
with an operator inserted. Finally, subscripts denote functional differentia- 
tion and * integration over the corresponding space-time argument. 

With the help of power counting we can identify the possible contributions 
to the right hand side of (|T0|) that have the correct tensor structure. It can 
be shown for all operators, with the exception of and Oj 2 \ that 
such contributions have to be proportional to Cp- The explicit calculation 

(2) (2) 

shows that this feature remains true for & and O) as well. Furthermore, 
there are no ^-dependent counter terms that do not vanish by the EOM. Our 
result (H) passes these consistency checks as well. 

4.2 Reparameterization invariance 

The HQET Lagrangian is a reparameterized form of the QCD Lagrangian, 
therefore the matching coefficients C of different orders in 1/rriQ are re- 
lated [TPl. For example the matching coefficient of the chromo-magnetic 



operators in 0(l/m 2 Q) can be derived from the coefficient in 0(l/rriQ) 

Cf ) = 2C { 2 1] - 1 . (11) 
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On the other hand, we know that the product of the matching coefficients 

and the renormalization constants CZ^ is finite and we can derive relations 
between the 0(a) matching coefficients, the tree-level matching coefficients 
and the anomalous dimensions: 

-» (a) -> (tree) 

C +C 7< 2 > = 0. (12) 

Since there are reparameterization invariance relations among the 0{a) match- 
ing coefficients, (|T2|) induces reparameterization invariance relations among 
the anomalous dimensions J7[. These relations are satisfied by our result (B|). 



5 Comparison with earlier calculations 

Two calculations of the renormalized HQET Lagrangian at order C(l/mg) 
have been circulated as preprints in the past || f|. Their results are not 
consistent with each other and our result differs from both. Therefore a brief 
discussion of the errors in these calculations seems to be in order: 

• After transforming the result of to our basis, it turns out that the 

~ (2) 

coefficients 7V of the local counter terms for the double insertions are 
incorrect. In particular, the — Ca entry in (jSc| ) is fixed by reparameter- 
ization invariance, which is therefore violated by the result in |§. The 
argument is unfortunately technically involved and will be presented 
elsewhere 0. 

• The operator basis used in Q is inconsistent. While these authors have 
used the EOM in their calculations, they do include an operator 

I (h v <7 a »T a h v ) D a F« u v» = - 9 -h v a^ ivD] h v (13) 

in their basis, which vanishes by this equation of motion. On the other 
hand, they have missed a spin-symmetric operator that is required as 
a counter term for the double insertions of the kinetic and chromo- 
magnetic operators. 



6 Renormalization group flow 

An analytical solution of the renormalization group equation using (f|) in the 
full basis seems to be impracticable. We can however restrict ourselves to 
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the basis of the operators O^T^ that do not vanish by the EOM: 
/ -\C A \ 



-(2) 
7phys 



ls~< 8 si 
"6^ ~~ 3°^ 



V 







-c A 








-\C A 








-c A ) 



This reduced renormalization group equation 

^»(,) + ^<? <2 '(,)^o 

with the initial (matching) conditions 



Cf\m Q ) 
Cf\m Q ) 



C^{m Q ) = C\ L> {m Q )C\ L \m Q ) = l 



= -1 



.(i) 



(!)/ 



c?K) = ^K)^K) = i 



can be solved analytically 



C?\») 



8C F 
~Ca~ 



a(m Q ) 



6/3« 5 



- 1 



.a) 





C A 




4/3« 


<x(™ Q )J 








a{fi) \ 









(14) 



(15) 

(16a) 
(16b) 
(16c) 
(16d) 
(16e) 

8C F 1 

c7~ 2 

(17a) 

(17b) 
(17c) 

(17d) 
(17e) 



Flavor threshold at which (3^ changes have been ignored in fllTD - It is 
straightforward to recover them by pasting solutions together at the thresh- 
olds. 

The solution of the renormalization group equation in the full basis can 
be obtained numerically for specific values of the gauge parameter £ and a 
specific gauge group (i.e. SU(3)). This will be done for the renormalization 
of the currents ||. 

7 Conclusions 

We have presented a new calculation of the renormalized HQET Lagrangian 
at order C(l/mg). Our result corrects previous calculations and obeys the 
Ward identities imposed by the BRST invariance of QCD and the reparam- 
eterization invariance of HQET. 

The new renormalized Lagrangian has been used in a calculation of the 
renormalized HQET currents at order C(l/mg). The results will be pub- 
lished in a sequel || to this note. A more detailed discussion of the consis- 
tency checks provided by BRST and reparameterization invariance will be 
presented elsewhere 0, together with technical details of the calculation of 
renormalized Lagrangian and currents. 
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